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Abstract 

The effect of non-zero c quark mass on b quark HQET Lagrangian, up to 1/mb 
level, is calculated at two loops. The results are expressed in terms of dilogarithmic 
functions of m c jm^. This calculation involves on-shell two-loop propagator-type 
diagrams with two different masses, m^ and m c . A general algorithm for reducing 
such Feynman integrals to the basis of two nontrivial and two trivial integrals is 
constructed. 
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1 Introduction 



On-shell two-loop calculations in any theory containing two massive fields with different 
masses and a massless field necessarily involve diagrams like Fig. |I[ Such diagrams appear, 
for example, in QED with e, \x and r, in QCD with b and c quarks, and in the electroweak 
theory. In general, such diagrams have a two-particle threshold^ at p 2 = M 2 , and a three- 
particle threshold at p 2 = (M + 2m) 2 , where p is the external momentum. There are 
two three-particle pseudothresholds, at p 2 = (M — 2m) 2 and p 2 = M 2 . Therefore, the 
on-shell condition p 2 = M 2 means that we are, at the same time, at the two-particle 
(pseudo)threshold and at the three-particle pseudothreshold. 



m 




Figure 1: Two- loop self-energy diagram. 



Combining the identical massless denominators, we can express these diagrams via 
scalar integrals (see Fig. ^) 



I(ni,n 2 ,n 3 ,ri4) = 



1 r d d k d d l 



(1) 



7T d J v^v^ 2 v^v2 4 

where 

V\ = M 2 — [p + k) 2 , V 2 = -k 2 , V 3 = m 2 - I 2 , V A = m 2 - (k - I) 2 , (2) 

and d = 4 — 2e is the space-time dimension (in the framework of dimensional regulariza- 
tion U). Evidently, 

I(ni, n 2 , n 4 , n 3 ) = I(ni, n 2 , n 3 , n A ) . 

Some integrals of this class have been already considered in the literature Here 
we propose a general algorithm, based on integration by parts Jl 1|| , which reduces any 



integral (P (with integer Uj) to the basis of two nontrivial integrals (known near four 
dimensions up to finite terms in e) and two trivial ones. 

In realistic calculations, we also get integrals ([!]) with scalar numerators. In some cases, 
these numerators cannot be directly expressed in terms of the denominators @). Such 
integrals are denoted as I(ni, n 2 , n^, n 3 ; no) (where n Q is the power of the corresponding 

1 It coincides with the corresponding pseudothreshold, since one of the particles involved is massless. 
2 The asymptotic expansion of similar on-shell integrals has been constructed in P,[lC|]. 
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p p + k p 

Figure 2: On-shell two-loop integral with two masses. 



numerator, see eqs. fll8|)-(|T9|)), and they are considered in Sec. 0. They can be reduced to 
the same basis as the integrals I(ni, ^2,^4,^3). 

The integration-by-parts technique [|11]] was first developed in the context of massless 
calculations. Later it was used in on-shell two- loop calculations with a single mass |3],f|, 
as well as in HQET jnj (see a short review [H)]). 



Recently it was shown |T6| (cf. also |17]]) that any two- loop propagator integral with 
generic masses and external momentum can be reduced to a finite set of basis integrals^. 
However, in special cases like the thresholds and pseudothresholds some general formulae 
of ref. ]l(J become degenerate. Therefore, these cases require a special examination. In 
particular, here we propose a reduction algorithm for the integrals ([!]) and implement it as 
a REDUCE |H||0| package. 



Integrals of this class appear in matching QCD to Heavy Quark Effective Theory 
(HQET, see |2l| for review and references). Specifically, let us consider the b quark HQET 
with rrib = M, keeping one more massive flavour, c, with m c = m. Coefficients in the 
HQET Lagrangian (as well as in the 1/M expansion of QCD operators) are derived by 
equating the on-shell matrix elements in QCD and HQET. In order to determine coeffi- 
cients in the Lagrangian up to the 1/M level, it is sufficient to equate scattering amplitudes 
of an on-shell heavy quark in an external gluon field: 




2M ^ ' AM 

= uMt a (v» + £j + cMz-^jiJ^gpj mo) + o . (3) 

Here p = Mv, q = p f — p, p 2 = p' 2 = M 2 ; form factors on the QCD side are F\(q 2 ) = 1 + 
0(q 2 ), G m (q 2 ) = fi g + 0(q 2 ); C m (/x) is the coefficient of the chromomagnetic interaction, the 
only nontrivial coefficient in the HQET Lagrangian to this order; Z _1 (/i)/igo is the matrix 
element of the renormalized HQET chromomagnetic operator (note that C m (/i)Z _1 (/x) does 
not depend on the normalization scale fi); the spinors are related by the Foldy-Wouthuysen 
transformation u(Mv + k) = (1 + #/(2M) + 0(k 2 /M 2 ))u v (k) (see |22| for a recent short 



3 The algorithm of jl6| was originally implemented in FORM and then p8[ in Mathematica. 
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review of the status of calculations of coefficients in the HQET Lagrangian). Two-loop 
anomalous dimension of the HQET chromomagnetic operator (or Z) was calculated in |f23| , 
24]. The chromomagnetic coefficient C m (fi ~ M) was calculated at two loops in |2l|, under 



the assumption that all other flavours (except the heavy flavour of HQET) are massless; 
higher orders in the large /3q limit were summed in |[25| . 

Here we calculate the effect of non-zero c quark mass on the chromomagnetic interaction 
in the b quark HQET. The HQET on-shell loop integrals without massive quark loops 
contain no scale and hence vanish. Therefore, only the integrals involving c quark loops 
are relevant. A simple method for calculating such integrals was proposed in || (and 
used in ||24j|). The QCD calculation involves the integrals ([!]). The non-zero m c effect on 
heavy- light bilinear quark currents in HQET was considered in ||. 



2 On-shell two-loop integrals with two masses 



Integrals ([!]) with n 3 < or n 4 < reduce to products of one-loop ones, and are propor- 
tional to T\ = 1(1,0, 1,0) with rational coefficients. Integrals with n\ < are two-loop 
vacuum "bubbles" , an explicit formula for them can be found in || ; they are proportional 
to T = 1(0, 0, 1, 1) with rational coefficients. These trivial basis integrals are 



T = (m 2 ) d - 2 T 2 {\ - d/2) , T x = (Mm)"r 2 (l - d/2) 



(4) 



Explicit general formulae for similar integrals with numerators can be found in ref. p6|P|. 

When rii, n 3 and n 4 are all positive, we employ the integration-by-parts technique |TT| . 
Applying the operator -^-(k—l) to the integrand of ([!]) and substituting p-l — > (p-k)(k-l)/k 2 
in the numerator, we obtain the recurrence relation 

4m 2 n 4 4 + I = \-2d+m + 2n 2 + 4n 4 - (n x + 2n 2 )(3" - 4~)2 + + ni(2~ - 3" + 4~)1 + ] / , (5) 

where, for example, 4 I(ni, n 2 , n 4 ) = 1(111,712,113,114 ± 1), etc. Due to the symmetry, 
we always may assume that n 4 > n 3 ; the index n 4 can be reduced down to n 4 = 1 by (|5|). 
After that, we are left with I(iii, n 2 , 1, 1), plus trivial integrals which can be expressed in 
terms of (§). 

Next we are going to reduce ii\. To this end, we consider recurrence relations obtained 
by applying ^ • k and • p to the integrand of ([[]) : 



d-ni-2n 2 -n A - ni2~l + - n 4 (2~ - 3")4 + 







(6) 

= 0. 

(7) 

We express 2 4 + I from @ and 1 4 + J from its 1 2 + shifted version, and substitute both 

into any. 



-2ni+2n 2 +n 4 +2ni(2-+2M 2 )l + -(2n 2 +n 4 )l"2 + +n 4 (l-3"2 + -l- + 2--3^)4 H 



d - 2m - 1 - (d - m - l)l-2 + + m(2 _ + 4M 2 )1" 



0. 



'Note that 7(— 1, ?i2, 713, 714) = 1(0, 712 — 1, «3, n^). 
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Now we express 4 + 7 from and substitute it into the 2 + shifted version of ([]); adding 
and substituting n 3 = n 4 = 1, we arrive at 



4(M 2 - m 2 )n x l + I(n x , ra 2 , 1, 1) 



-3d + 3n x + 2rt 2 + 5 + (d - m - 1)1"2 + 
4m 2 (3~4+ + d - rii - 2n 2 - 3)2+1 J(m, n 2 , 1, 1) .(9) 



This relation allows one to reduce n\ down to n% — 1 (note that the term with 3~4+ is 
trivial here). Thus we are left with 7(1, n 2 , 1, 1). 

Finally, substituting 1 + J from @ and its shifted version 2 l + 7 into (§), we obtain at 
ni = n 3 = ri4 = 1 

(3d - 2n 2 - 6)2" + 4(M 2 + m 2 )(2d - 2n 2 - 5) + 16M 2 m 2 (d 

= (d - 2)l~(4m 2 2+ + 1) - 4m 2 3"4+(4M 2 2+ + 



2n 9 



4)2+ 7(1, n 2 , 1,1) 



J l,n 2 ,l,l . 



(10) 



This relation allows one to lower or raise n 2 (note that the integrals on the right-hand side 
of eq. (|T0| ) are trivial). Therefore, all integrals (|l|) can be expressed via 



J = 7(1,0,1,1), 



h 



7(1,1,1,1) 



'ir 



and the trivial integrals To,i (H) ; exactly in d dimensions. The basis integrals (|TT|) in d 
dimensions can be expressed via hypergeo metric functions 37^, see Appendix A. 

Expressions for 7^1 expanded in e up to the finite terms can be, with some efforts, 
extracted from any two independent integrals out of those calculated in 0-0 • Some details 
of this procedure are discussed in the Appendix B. Introducing a dimensionless variable 



r = m/M 

and dilogarithmic functions^ 
L + = - Li 2 (-r) + 1 log 2 r- logr log(l+r) - \ix 2 = Li 2 (-r _1 ) + log r" 1 log(l 



(12) 
(13) 



L_=Li 2 (l -r) + |log 2 r + 



1^2 



Li 2 (l 



Li 2 (r) + ilog 2 



logr log(l — r) + \n 2 



=Li 2 (r 1 )+logr 1 log(l — r l ) (r > 1 
the results for the integrals 7 01 can be presented as 
7 



T 2 (l + e) 



7i 



M 



2-4e 



— m 



2-4e 



M' 



-de 



■ 1 5 

1 3 

H 2 logr + 6 



5 



„2\2 



1^2 
77 7T 



(r < 1) 
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(14) 



(L + + L_) -21og 2 r + — 



+ 0(e) 



-^ + - + 2!1 +,-)-'L + + 2(l - -rfL. 
2e z 2e 



nl 2 19 

21og 2 r + y 



+ 0(e) 



(15) 
(16) 



5 We adopt the notation L± which has been used in [g]45| (cf., e.g., eq. (A3) of ||). In Q, the functions 
i_ and L + were called i?i and i?2, respectively. In || the functions T + (r) = — 2L_ + ^ log 2 r + ^ir 2 and 
ilog 2 ; 



T-(r) = — 2L_| 



7r were used. They have a nice property T (r ) = — T (r), and they are 



analytic continuations of each other atr<-> — r. 
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Eq. fll5|) corresponds to a special case of more general results for the sunset-like diagrams 
presented in ||. Eq. (|^) is equivalent to the result presented in eq. (29) of ref. ||. In 
Appendix B, we discuss different ways to derive it. Note that 

L + + L_ = | Li 2 (l - r 2 ) + log 2 r + ^vr 2 = -§ Li 2 (l - r" 2 ) + ^vr 2 . (17) 

When dealing with diagrams like Fig. |IJ with numerators and reducing them to scalar 
integrals, one cannot directly express I -p via Pi. ..4. Therefore, integrals similar to (fj) but 
containing some extra numerators^ appear. Let us define 

~ 1 f d d kd d l M no 

I{n 1 ,n 2 ,n 3 ,n4;n ) ■ 



(18) 

p 2 =M 2 



7T d J V^V^V^V^ 

As compared to (|I|), it contains an extra factor M n ° in the numerator, with no > and 

M= (21 - k) -p . (19) 
Changing the integration momenta I <-» k — I, we obtain 

I(ni,n 2 ,n4,n 3 ;n Q ) = (-l) n °I(ni, n 2 , n 3 , n 4 ; n ) . 

In particular, when n 3 = the integrals with odd uq vanish. Integrals with < or 
ni < are trivial, as before. 

In principle, one could deal with (|1|) by generalizing (to higher uq) the substitutions 
like p ■ I — > (p ■ k)(k ■ l)/k 2 (cf. eq. (^)). However, the higher no is, the more cumbersome 
these substitutions are. An alternative, recursive way is to use the integration-by-parts 
procedure^. Applying J| • p to the integrand, we obtain the recurrence relation 

n 4 4 + + J = f-7ii+7i 2 - M 2 ri CT +rix(2~ + 2M 2 )1 + - n 2 l"2 + - fn 4 (l~ - 2~)4 + l J. (20) 

We may assume n 4 > n 3 ; if n 4 > 1, then both n and n 4 can be reduced by (^). Thus we 
are left with already known integrals and I{n\, n 2 , 1, 1; no) with even no- Applying Mr ■ I to 
the integrand, we obtain at n 3 = n 4 = 1 

n 1 l + + I(n 1 ,n 2 , 1, l;n ) 
■\n (l- - 2 r)0- + 2nx3-l + + 2n 2 3-2+ + (2" - 3~ +2m 2 )4+j I(n 1} n 2 , 1, 1; n ) .(21) 



If ni > 1, then both no and ni can be reduced by (|21|). Now, we are left with known 
integrals and 1(1, n 2 , 1, 1; n ) with even n . Finally, applying • Z to the integrand, we 
obtain at ni = n 3 = n 4 = 1 

(d + n - 3)/(l, n 2 , 1, 1; n ) = -(2" - 3" + 4m 2 )4 + /(l, n 2 , 1, 1; n ) . (22) 

The right-hand side contains integrals with the same no but n 4 > 1; in them, no can be 
again reduced by (|20"1). In this way, these more general integrals (|18D are also reduced to 
the basis ones, J 0j i and T 0jl . 

6 These numerators are not really irreducible, like e.g. in the three-point functions j27j], since the corre- 
sponding integrals can be reduced algebraically. 

7 A similar algorithm has been constructed in fwh. 
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3 Effect of m c on b quark chromomagnetic interaction 



First, we calculate contributions of a flavour with mass m (e.g., the c quark) to the mass 
and wave function renormalization constants of a quark with mass M (e.g., the b quark) 
in the on-shell scheme (Fig. |l|) in d dimensions. The results are[| 



AZ 



M 



TfCf 



9 4 o 



(4vr) d 



(d-2) 2 (d-5) 
(d-3)(d-Q)M i \2r 2 

9o d 



d 2 -9d+lQ T 4(d-5)(l+r< 



AZr 



T F Cp 



(An) d 1 - r 2 



(d-6)M ; 



-Io + 



d-6 



-h 

(23) 
(24) 



i=i 



where 



(d - 2) 2 T 



(d-3)(d-5)(d-Q)M 2 m 2 



Jo 



{d-2) 2 T 1 



(d- l)(d - 3)(d - 6)M 4 



J\ 



h 



(d- l)(d-6)M 2 ' (d-l)(d-6)' 
The results for the coefficients a, in AZq are 

oi = (2d 4 r 4 - d 4 r 2 - 39d 3 r 4 + 21d 3 r 2 - 2d 3 + 272d 2 r A - 159dV + 30d 2 

- 789dr A + 509dr 2 - lUd + 770r 4 - 566r 2 + 216)/[4(d - 2)(d - 7)r 2 ] , 
a 2 = (d 2 + 5dr 2 - 12d - 25r 2 + 31)/2 , 

a 3 = {2d 3 r 2 + d 3 - \7d 2 r 2 - I2d 2 + 27dr 2 + A7d + 8r 2 - 56)/2 , 
a 4 = 2(d 3 r 2 - 12d 2 r 2 + d 2 - 5dr 4 + 45dr 2 -6d + 25r 4 - 54r 2 + 5) . 

After expansion in e up to finite terms, these formulae reproduce the results of ||,£§. 
In the case m = 0, 



T = Ti = , h 



3d -8 
A(2d-7) 



M 2 I . 



(25) 



Note that J at m = is called -M 2d ~ 6 h in ||]. AZq contains a term Ji/r 2 ~ T /m 4 
which does not vanish in the limit m — > 0. This means that there is no smooth m — > 
limit in AZq [ffl. Omitting the T 0j i terms and then setting m = 0, we reproduce the results 
of for the massless quark contributions, exactly in d dimensions. 
In the case m — > M, 



Ti = T 
h ~- 



d-2 



2 

3d -8 



1-r 2 ) 



A(d-A)M 2 
+ 0((l-r 2 ) 2 



1 



+ 0((l-r 2 ) 2 ), 
■ 2 l 3(d 



8(d-3)(d-4)M 4 



r 



2d + 3, 
1 + — 6~ (l-O 



(26) 



'Id QCD, T F = l,C F = (N 2 - 1)/{2N), C A = N, where N = 3 is the number of colours. 
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21 



M 2d ~ 4 ! 2 and 



-M 2d - 6 I 2 , 



respectively. 



Note that To and Jo at m = M are called in 

Using this, we reproduce the results of ]3|,f| for the contributions of the quark with mass 
M, exactly in d dimensions. 

Now we calculate the contribution A/x g of a flavour with mass m (say, c) to the chromo- 
magnetic moment of a quark with mass M (say, b) (Fig. |3|), taking into account the wave 
function renormalization fl24|). We use the background field method | 28fl . Discontinuity 
in the limit m —>■ in the Cp term has cancelled with Zq, but in the Ca term it is still 
here. We reproduce the m = result (by omitting T 0j i terms) and the m = M result from 
Appendix in [p4j| , exactly in d dimensions^. 




(a) 



(b) 



(<0 



Figure 3: Diagrams for the chromomagnetic moment: (b) implies also the mirror symmetric 
diagram, and (c) has two orientations of the quark loop. 



Combining this result with the corresponding HQET term A/i g0 W% i we obtain 

4 4 



A/i g — A/i g o = Tj 



9l 



(47r) d 1 



Cf X! a FiJi + Ca X! a Ai Ji 



i=l 



i=l 



(27) 



with the coefficients 



a F1 
a F2 

CtF3 

a Pi 

dAl 

aA2 
aAi 



{d 3 r 2 - lldV - 3d 2 + 33dr 2 + 30d - 19r 2 - 71) /2 , 
-(rf 3 r 2 - 2d 3 - YId 2 r 2 + lid 2 + 89rfr 2 - U6d - 133r 2 + 151) , 
d 4 r 2 + 2d 4 - 7d 3 r 2 - 31d 3 - 7d 2 r 2 + 168d 2 + 105dr 2 - 375c/ - 152r 2 
4(d 4 r 2 + d 3 r 4 - 15d 3 r 2 - 17d 2 r 4 + 85d 2 r 2 - d 2 + 89dr 4 - 22\dr 2 + 6d 
- 133r 4 + 210r 2 - 5) , 

-(d 3 - d 2 r 2 - 13d 2 + Wdr 2 + 52d - 26r 2 - 64) /8 , 

(2d 3 r 2 - 3d 3 - 29d 2 r 2 + 39c/ 2 + 129c/r 2 - 156c/ - 182r 2 + 200) /4 , 

-{3d - 8)(c/ 3 + d 2 r 2 - 11c/ 2 - Wdr 2 + 38c/ + 26r 2 - 44)/4 , 

-(c/ 4 r 2 + 2c/ 3 r 4 - 18c/ 3 r 2 + 2c/ 3 - 29c/ 2 r 4 + 114c/ 2 r 2 - 18c/ 2 

+ 129c/r 4 - 299c/r 2 + 44c/ - 182r 4 + 282r 2 - 28) . 



296 



Discontinuity in the limit m — > in the Ca term has cancelled. The difference of this 
contribution and that of a massless flavour is finite as e — > 0; it represents the amount by 



'There is a typo in the Appendix of p4|: Jo, 1,2 should read Ji,2,: 
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which the b quark chromomagnetic interaction coefficient changes due to a non-zero mass 
of c quark:[|] 



AC m — T> 



An) 



8C F A F + ^C A A A 



(2f 



Ap = -r(l + r)(l - r - 4r 2 )L + + r(l - r)(l + r - 4r 2 )L_ + 6r 2 (logr + | 
A A = -(1 + r)(2 + 4r - r 2 )L+ - (1 -r)(2 - 4r - r 2 )L_ + 2 log 2 r + ±7r 2 + 2r 2 (logr + 1) . 

It vanishes as r — > 0, whereas at r = 1 it reproduces the result of ref. (namely, the one 
with heavy quark loop). 

In the QED case (Tp = 1, CV = 1, CU = 0, a s — > a), the formula (£8|) is closely related 
to the contribution of a lepton with mass m (e. g., electron or r) to the magnetic moment 
of a lepton with mass M (e.g., muon). Adding back the massless contribution (which was 
subtracted in fl28|)) and re-expressing the MS coupling in the one- loop contribution via the 
on-shell coupling a(M) = a (l — |^logrJ, we find that this contribution is given by ( p8|) 

with Ap — > Ap — | logr — ||. It was calculated in [0]; our result agrees with it (and has a 
simpler form). 

Expanding (p8|) at r 1 and r > 1, we obtain 



Ap = Wr + 2(4 log r + 3)r 2 - fvrV + 2(2 log 2 r - ^ log r + ±7r 2 + f ) r 4 



2 

oo 



+ E I 2 9F (2 n ) log r + ^ ) r- = £ ( -2 9P (-2„) log, + iMzM 

n=3 \ / n=0 \ 



r -2n 



CXJ / 

A A = 3vr 2 r + 3 (log 2 r + 4 logr + |tt 2 - 3) r 2 - ±vr 2 r 3 + ^ 2# A (2n) logr + 



n=2 



2 log 2 r + f logr + |vr 2 + f + £ f-2^(-2n) logr + 2n) 



n=l 



dn 



r 



dgu(2n) 
dn 

-2n 



r 2n 



1 5 4 , , 2 6 3 1 

9f{x) = + , g A {x) = — + + . 

x — 1 x — 3 x — 4 x x — 1 x — 2 x — 3 

Note that loop effects of a very heavy flavour (say, t) do not decouple in C m , because it 
is not directly measurable. Adding (— |logr — |§) to the expansions of Ap, we obtain 
the expansions of the contribution of a lepton with mass m to the magnetic moment of a 
lepton with mass M in QED; they agree with the expansions of obtained in It is 
also easy to obtain as many terms as needed in the expansions as r — * 1 using the formulae 
for L_ via 1 — r or 1 — r _1 (see eq. flT4])) and for L + + L_ via 1 — r 2 or 1 — r -2 (see eq. QT7D). 



4 Conclusions 

We have presented an algorithm which reduces any scalar integral I(nx,ri2, n^, n^, no) (given 
by eq. (|I8D) to two trivial basis integrals (|J) and two nontrivial ones, whose e-expansions up 

10 In the notations of |5), A F = A 2 - 4A 3 . 
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to the finite terms are given by eqs. ( |T5D and flTo]). The algorithm has been implemented as 
a REDUCE package. It has been tested by evaluating about 27000 recurrence relations for 
specific values of rij, including those relations which were not directly used for construction 
of the algorithm. This package allows one to calculate, completely automatically, the on- 
shell two-loop self-energy diagrams (Fig. |1]) in any theory with two massive particles having 
different masses^. 

We calculated the effect of non-zero m c on the b quark chromomagnetic interaction 
coefficient fl28|). Combining it with the result for m c = [^], we have numerically 



C m (m b ) ~ 1 + lj W m& ) + (14,1439 + A m {m c /m b )) (-Y , A m (0.3±0.03) = 0.98±0.07. 

6 7T V 7T / 

(29) 

The effect is moderate and positive. This is the only nontrivial coefficient in the HQET 
Lagrangian up to l/m b level. 



Due to reparametrization invariance |30[ , it also determines [ 3TJ] the spin-orbit interac 



tion coefficient — the most important spin-symmetry breaking coefficient at 1/m^: 

CM = 2C m (ji) - 1 . (30) 

This relation has a very simple physical interpretation.^ The spin-orbit interaction is the 
interaction of a moving (nonrelativistic) heavy quark with chromoelectric field. In the 
quark rest frame, the field acquires a chromomagnetic component. The quark chromomag- 
netic moment interacts with it, producing the term 2C m . There is also Thomas precession, 
which compensates half of the chromomagnetic contribution at the tree level. This is a 
purely kinematic effect due to Lorentz transformations, and hence it gets no corrections. 

The influence of m c on other l/m b terms in the b quark HQET Lagrangian, as well as 
on l/m& expansions of various QCD currents, can also be calculated using the presented 
method. 
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Appendix A. Integrals in d dimensions 

Here we present the most general results for the integrals ([!]). To indicate that the results 
of this appendix are valid for arbitrary powers of propagators, not only the integer ones, 
we substitute rij — ► v^. 

11 The package may be obtained from http://www.inp. nsk.su/^grozin/. 
12 AGG thanks I. B. Khriplovich for this interpretation. 
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Using Mellin-Barnes contour integral representation for the massive denominator T> 1 1/1 
(for details, see in [0]), we get 



^(^1,^2,^3,^4) 



r (z^ r(z^ r(z/ 4 

— 1CX) 

r(i/ 3 + s) I>4 + s) r (1/3 + z/ 4 - f + s) T(d -v x - 2v 2 + 2s) 



x 



T(z/ 3 + z/ 4 + 2s) T(d — V1 — V2 + S) 

where the integration contour separates right and left series of poles of V functions occurring 
in the numerator of the integrand. 

If we close the contour to the right, we should sum over two series of polesQ. The 
general result can be presented as 

,2\d— v\— f 2 -f3-^4 



^2, ^3, ^4) 

r(f-Z/i-l/ 2 ) r(z/ 1 + I/ 2 + Z/ 3 + Z/ 4 -c/) r(z/! + Z/2 + Z/3- f) r(z/!+z/ 2 + 



X 



+ 



i> 3 ) r(z/ 4 ) r (f) r(2z/i + 2z/ 2 + ^ + z/ 4 - d) 

^1 + ^2 + ^3-f, ^1 + ^2 + ^4-f, Z/1+V2+V3 + I/4 — d, 5^1, + 1 

^l + ^ + K^ + ^-f, Z/l + Z/ 2 + |(z/ 3 + ^4+l)-|, l/l+i/ 2 -f + l 

r (1/1 + r (i/3 + 1/4 - f) r(d - 1/1 - 2v 2 ) 



x 5 F 4 



' 2 



M 2 ' 
m 2 



'M 2 \ 2- Ul ~ V2 



X 5F4 



r(i/i) r(z/ 3 + 1/4) r(d - ^ - 1/ 2 ) 

^3, ^4, ^3 + ^4 - 5, f (d -v\- 2z/ 2 ), §(d - V\ - 2v 2 + 1) 
§(z/ 3 + 1/4), |(^3 + 1/4 + 1), d - r/ x - z/ 2) § - 1/1 - v 2 + 1 



M 2 ' 
m 2 



,(31) 



where 



/ ai, . 


. ,a P 


•) 


V c i, • 


•,cq 





E 

i=o 



a P 



3 ' ' ' 



{CQ)j f- 



r(a + j)/T(a) is the Pochhammer sym- 



is the generalized hypergeometric function; 
bol. 

For trivial cases 2/1 = and 1/3 = (or 1/4 = 0), the expression (|3T| ) reproduces known 
results in terms of T functions. When v 2 = (the "sunset" configuration), the same 
result (|31~D can be reproduced in this limit (i.e., p 2 = M 2 and two other masses are equal) 
by taking two sums in the hypergeometric series of three variables presented in ref. |34| . 

One can see that for integer (non-negative) values of Vi both 5F4 functions in (^) reduce, 
for arbitrary d, to a finite sum of 3-F 2 functions of the same argument. For example, we 
get for our basis integrals 7 

h 



7(1,0,1,1) and h 

2\l-2e r y 



7(1,1,1,1) 



[m 



e 2 (l-e) 



3 r 2 



1 I 

■ L ' 2' 



-1 + 2e 



e, I + e 



M 2 



v 1 2 

13 If we close the contour to the left, we get a more cumbersome result in terms of the hypergeometric 
series of m 2 /M 2 , corresponding to the analytic continuation of (pl|). 
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exactly in d = 4 — 2e dimensions^. Expansion of Jo and I\ in e up to finite terms is given 
in eqs. (]T§)-(]T§). 



Appendix B. Various ways to calculate basis integrals 

Using a procedure similar to M, the result ([Tj]) can be derived as follows. Note that m i£ Ii 
and m 4£ /(2,0, 1, 1) are dimensionless and depend only on r = m/M. Differentiating the 
definition (|IJ) in M and taking into account that p = Mv (v 2 = 1), we see that I± should 
obey the differential equation^ 

o d r 



dr 



Therefore, 



r- 1 m Ae h\ = -m 4e /(2, 0,1,1) 
l: r / d P 



h=r j-fm 4e 7(2, 0,1,1) 



lr— >p 



(32) 



(33) 



The expression for 1(2, 0, 1, 1) can be taken, e.g., from ref 
1(2,0,1,1) „ 2E 



r 2 (l + e) 



(M 



^ + ^-^-21og 2 r + 2(l-r 2 )(L + + L_; 



+ 0(e) . (34) 



Substituting eq. ([34]) into ([33]) we get (|16|). 

An alternative derivation follows the method used in 0. The integrals ([I]) can be 
presented as 



I(ni,n 2 ,n 3 ,n 4 ) 



i r U 12 (k 2 )d d k 



71 



d/2 



with n n3n4 (fc 2 



p 2 =M 2 



V n iv n 2 

For instance, when we consider integrals with = 1 and ?i4 = 2 we get 

A 2 -l 



d d l 



n d/2 J pnsp«4 • 



n 12 (fc 2 



4m 2 A 



log^i + 0(e), with A 
A — 1 



1 



4m^ 



(35) 



(36) 



14 If M — m, one arrives at a representation in terms of 3-F2 functions of unit argument |E 
15 This equation can also be derived using the method of [p3|. 
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Moreover, we know that Ui 2 (k 2 ) may be represented via dispersion integral, IIi 2 (&; 2 ) = 
/ dsp(s)/(s — k 2 ), the explicit form of the spectral density p(s) being not relevant for our 
discussion. 

As the first example, let us consider the integral 7(2, 0, 1, 2), which is convergent. Sub- 
stituting the spectral representation for II12, combining the denominators T>\ and (s — k 2 ) 
by the Feynman formula and, then, calculating the integrals over k and s, we arrive at 

7(2,0,1,2)= fdx-^- U l2 {k 2 )\ Ml ^ + 0{e). (37) 

J 1 X i—x 

If we consider, say, 7(3, 0, 1, 2), we get in the integrand (dU 12 (k 2 ) /dk 2 )\ k2= _ M2x 2/^ x y The 
second example is 7(2, 1, 1,2). In this case, the on-shell singularity can be separated via 
U 12 (k 2 ) = [n 12 (fc 2 )-n 12 (0)]+n 12 (0). Note that [n 12 (P)-n 12 (0)]/A; 2 = fdsp(s)/[s(.^k 2 )]. 
In this way, we obtain 



^(2, 1,1,2) 
T 2 (l + e) 



(Mm 



,-2-2e 



is 



dx 



x 



rii 2 (£; 2 )-IT 12 (0) 



x 



k' 2 



k 2._ 



M 2 x 2 



+ 0(e) . (38) 



Integrals with n 2 > 1 require more subtractions, and those with n\ > 2 involve derivatives. 
After inserting the expression (|36|), the integrals over x are of the following form: 



A(x) . X(x) + 1 _ . , 
X(x) \{x) — 1 



dx . 



(39) 



where X(x) = X 



k 2 =-M 2 x 2 /(l-x) 



(cf. eq. (BB)), whilst A(x) and B(x) are rational functions. 



Introducing the new variable y (the limits of the y integration are also from to 1), 
x(X(x) — 1) _ l1 r(l — y 2 ) , 1 dx r(l — y 2 ^ 



y 



so that x 



and 



2r y + r X(x) dy 

we see that all occurring structures become rational, since 

1 + 2ry + y 2 X(x) + 1 1 + ry 



(y + 



(40) 



X(x) 



y 1 



X(x) - 1 y(r + y) 



(41) 



Therefore, the integral ( p9|) can be expressed in terms of dilogarithms. In this way, 
from (|37T), (|38f) we obtain 



7(2,0,1,2) 



L+ + L_ 

M 2 



0(e), 



7(2,1,1,2) 
r 2 (l + e) 



1 



4M 2 m 2+4e 



-^+r(L + -L_) + 2 



■0(e). 
(42) 

This is enough to reproduce (|i~5D , ([161) . 

Similar procedure can be also applied to 7(ni, n 2 , n 3 , n 4 ) with other values of n i; as soon 
as a suitable subtraction can be performed. Apart from subtracting il n3n4 (0) (and higher 



terms of the Taylor expansion in k ), subtractions of II 



n 3 n 4 \ m . 



=0 may be also considered. 



13 



Some further results of interest are 
1(1,0,2,2) " L+ + L - 
7(2,0,2,2) 



L + -L. 



M 2 



AM 2 m 2 



r(L H 



0(e), 

2r 2 log r 
1 — r 2 



0(e). 



In particular, if we take 7(2, 0, 1, 2) and 7(1, 0, 2, 2) as independent finite integrals, we get 
only "pure" L + ± L_ combinations, without logarithms, etc. 

The results ( |T5| ) and ( |T6| ) can also be checked by using the formulae of Appendix A 
in ||. The quantities were calculated there, which are certain combinations of the 

integrals 



(M ) 



2\2— n\— ii2 



2 J(n 1 ,n 2 ) 



7T 



d/2 



P(k 2 ,m)-P(k 2 ,^ dk 



Am 2 



7(ni,n 2 , 1, 1) 



lm=0 



- 7(m, n 2 , 1, 1) + 27(m, n 2 + 1, 1, 0) + — — I(n u n 2 + 1, 1, 1 



where 



P(k 2 ,m) 



Am 2 



(d-2)k 2 



n n (fc 2 



m 



-2c 



rd + e) 



--+ 1 + 



d-2 



k 2 
2m 2 \ 



(43) 



n 



10 



k 2 



A log 



A + l 
A - 1 



+ 0(e) 



(44) 



is proportional to the fermion-loop contribution to the gluon polarization operator^]. When 
an integral contains neither ultraviolet nor infrared (on-shell) singularities, it can be cal- 
culated in four dimensions. Averaging over the directions of k in four-dimensional 
Euclidean space, we obtain A^/V^ 1 , where A = A_i = 1, A\ — x, A 2 = x 2 /(2 — x), and 
k 2 = —M 2 x 2 /(1 — x). Therefore, convergent integrals are given by 



J(ni,n 2 ) 



P(k 2 ,m) - P(k 2 ,0) 



fe 2 = 



,2^,2 dx . 



(45) 



with $ nin2 (x) = A ni (x)x 3 ~ 2ni ~ 2n2 (l - x) ni+n2 - 3 (2 - x) . 
In terms of these integrals, the quantities A 1 2) 3 used in can be expressed as 

Ax = | [27(1, l) + 7(2, 0)] , A 2 = | [7(1, l)-7(2, 0)] , A 3 = | [7(0, 1)-7(1, 0)-7(2, 0)] . 

(46) 

This coincides with the integral representations ||. Explicit expressions for A, given in 
eq. (A2) of [§[] (which can be obtained using the substitution (|40D provide three independent 
checks on our basis integrals. Note that the integral calculated in || is A = |(Ax + A 3 ). 

16 In refs. ||§ the notation n(-m 2 /k 2 ) = [P(k 2 , m) - P(k 2 , 0)] rf=4 has been used. To get the d- 
dimensional unrenormalized contribution to the gluon polarization operator, eq. ([3]) should be multiplied 
by-2T F g 2 {^y d l 2 {d-2)/{d-l). 
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The integrals^ 7 ] 

(M*^-<*7(r, n)~ ' f P(P ' m) - P( °' m) d d h - ^ - 1} T(n n 2 0) 
(M ) J(m, n 2 ) = --^j ^-n 2 d k - ^Z^jltni, n 2 , 2, 0) 

Am? 

- J(m, ra 2 , 1, 1) + 2/(72!, n 2 + 1, 1, 0) + ^T^K, "2 + 1, 1, 1) (47) 

have better infrared convergence but worse ultraviolet one. We define Aj ^,3 via 

J(2,1)=A 1; J(2,0)- J(1,1)=A 2 , (48) 
J(0, 1) - J(l, 0) - J(l, 1) = A 3 + 2^7^(0, 2) . (49) 

The combination of integrals on the left-hand side of fl49|) is rather interesting. As k — > 00, 
the combination of their denominators behaves as 

1 1 1 £ ( 1 

+ 77 



M 2 £> 2 M 2 Pi V{D 2 (2- £ )fc 4 VA; 6 

where (on the right-hand side) terms vanishing when averaged over the directions of k are 
omitted. The term 0(l/k e ) gives an ultraviolet-convergent integral. The term 0{e/k^) 
might have been omitted, if an arbitrarily large ultraviolet cutoff were considered. With- 
out the cutoff, it gives a purely ultraviolet contribution (similar to the axial anomaly) 
eJ(0, 2)/(2 — e). The term A 3 can be calculated at d — 4 using the convergent integral 
representation with $ i {%) — $10 {%) — ®u(%)- Using eq. (10), we obtained 

Ai = ±(l + r)(4-r + r 2 )L+ + ±(l-r)(4 + r + r 2 )L_ + (§ + ±r 2 ) logr + ^ + ±r 2 , 
A 2 = 3r(l - r 2 )(L+ - L_) + 2(2 - 3r 2 ) logr + f - 6r 2 , 

A 3 = 3r(l + r) 2 (l - 2r)L + - 3r(l - r) 2 (l + 2r)L_ - 3(1 - 4r 2 ) logr + f - 15r 2 

and checked that these expressions agree with <i-dimensional results obtained by our pro- 
gram as d — > 4. The integral calculated in § is A = ^(4Ai + 3A 2 - 2A 3 ). Eqs. (\T§ 
and flTBj ) can be reconstructed from any two out of four independent integrals from | 
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